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ABSTRACT 



A three-component balance system for use in measuring 
drag and lift forces as well as pitching moments acting on 
a wind tunnel model is presently available in the Department 
of Aeronautics. The original manual counterweight readout 
procedure was considered obsolete from the standpoint of 
data acquisition; therefore, a revision was incorporated 
using strain gauge load cells. The design and installation 
of the electrical readout system was accomplished. Cali- 
bration procedures were devised and performed in order to 
develop a loading output matrix for use with a rapid data 
acquisition system and compatibility with an analog-digital 
computer program. 
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I. INTRODUCTION 



The Aeronautical Engineering Department of the Naval 
Postgraduate School possesses two subsonic wind tunnels for 
use in the Department curricula and for student research. 

The two-stage drive tunnel has a maximum velocity of two 
hundred miles per hour and a test section that measures 
three and one-half by five feet. The single-stage drive 
tunnel also has a maximum velocity of two hundred miles 
per hour and has a test section size of thirty-two by forty- 
five inches. 

Any requirements for force measurement in wind tunnels 
must be met through the use of one AEROLAB 543 wind tunnel 
balance. Although this balance is capable of measuring 
both longitudinal and lateral forces with a high degree of 
accuracy, it is a beam type balance that requires mechanical 
manipulation in order to read the loads on the model. Aero- 
dynamic coefficients must be calculated by using additional 
information acquired from other sources. For these reasons, 
final results are not readily available for use in data in- 
terpretation during the course of the experiment. 

It would be highly beneficial to any person conducting 
a wind tunnel test to have at his disposal immediate exper- 
imental results. This could be realized through the use of 
a rapid data acquisition system which would analyze balance 
loading, tunnel dynamic pressure and model size parameters. 
The use of a digital computer program which would require 
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reading-in of the above variables would be the primary 
method for rapid calculation of test results. Access to 
the program through the use of a terminal at the testing 
station would offer great flexibility. Many of the com- 
ponents for this type of operation are available in the 
Department of Aeronautics, such as those provided by Collier 
(Ref. 1) and Weinzapfel (Ref. 5). The integration of the 
complete data system would result in a routine pi eduction 
procedure . 

The Aeronautical Engineering Department acquired a TASK 
MK I three component beam balance in October, 1958; the 
balance, for various reasons, has not been used for its in- 
tended purpose for a considerable length of time. This 
thesis serves primarily as a report on a design modification 
and testing of this balance as a proposed component of a data 
acquisition system such as described above. A second purpose 
of the thesis is to serve as an operational manual on the 
modified balance f< 



Task MK I Modi- 



FIGURE 1 



fied Balance 




10 





II. SYSTEM DESIGN 



A. FUNDAMENTAL OBJECTIVES 

The main objective of the project was the adaptation 
of the Task Corporation MK I three component beam balance 
sd that it would provide consistent electrical output sig- 
nals in order to form the principal element of a rapid data 

% 

acquisition system for use in wind tunnel testing. The out- 
put was to be linear over the desired load ranges of the 
balance. These load ranges were chosen as five hundred 
pounds in lift, seventy-five pounds of drag and seventy- 
five foot pounds of moment; the desired ranges included 
both positive and negative values and were derived by con- 
sidering a dynamic pressure of sixty pounds per square foot 
and a model wing area of three square feet in conjunction 
with a maximum lift coefficient of three, a maximum drag 
coefficient of four-tenths and a maximum pitching moment 
coefficient of five-tenths. 

The electrical adaptation was required to have the 
capability to provide zero references for the longitudinal 
forces and angle of attack with a model in place on the 
balance. Provisions had to be made to ensure consistent 
output values under possible variances of power supply and 
other random disturbances. 

B. AN INTERIM PROPOSAL 

A proposal that was believed to meet the above stated 

i 
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objectives was formulated previously. The design utilized 
reluctance gauge load cells to measure the forces in the 
internal mechanism of the balance. The deflections of the 
load cell rings were sensed by the reluctance gauges, which 
in turn transmitted electrical signals whose strength was 
proportional to the forces applied to the load cells. 

The design eliminated the counterweight force measurement 
system of the original balance and appeared to be a reason- 
able solution. Once the design was fabricated, however, it 
was determined that the reluctance gauges were too sensitive 
to random disturbances and under applied loads produced un- 
acceptable nonlinearities in their output signals. After a 
length of time spent in investigating these problems, the 
proposal was determined to be impractical. However, the 
reluctance gauge load cell design did verify that the use 
of electrical damping in a voltage signal that oscillated 
about an average value was feasible. In a mechanical bal- 
ance as an analogy, viscous damping is provided through 
the use of adjustable dashpots in order to remove random 
load disturbances. 

C. A FINAL PROPOSAL 

Although the interim proposal was unusable because of 
the problems associated with the reluctance gauges, the con- 
cept of measuring forces applied to the balance through the 
use of load cells attached to its internal mechanism offered 
a prospective solution. The load cells were therefore 
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redesigned to use electrical strain gauges to measure load 
cell deflections. The new load cells were then fabricated 
and installed so as to replace the reluctance gauge load 
cells . 

The resistance elements of the load cells formed a 
portion of a Wheatstone bridge circuit with the differential 
voltage output being handled by contemporary operational 
amplifier circuitry. Each bridge was designed to be pow- 
ered by a five volt direct current supply, which was 
closely controlled by solid state regulators. The accuracy 
of these voltage regulators eliminated the need for separate 
adjustable power supplies for each bridge circuit. 

A control box (figure 3) was fabricated which enclosed 
the power and signal conditioning circuits for the load 
cells. The face of the box contained four bridge balancing 
potentiometers for the three load cells and the angle of 
attack transducer attached to the angle of attack drive 
motor. The control box was powered by 110 VAC and has 
power/signal cables with quick disconnect fittings for at- 
tachment to the three load cell bridges and the angle of 
attack transducer. 

The control box serves as a distribution point so as 
to route the output signals of the four bridge circuits to 
the multiplexing unit of a Digitec data acquisition and re- 
cording system (See figure 21) . This arrangement allows 
rapid recording of four channels of output information 
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whose values are a measure of lift, drag and moment loads 
applied to the balance as well as of the angle of attack of 
a model . 

In the final application, an additional signal pro- 
cessing unit would be employed which includes provision 
for electrical damping, scale adjust and signal monitoring. 
This unit was originally designed for use with the reluc- 
tance gauge load cell system and is intended as an inter- 
face between the load cell signal conditioning amplifiers 
and the output paper tape punch device that is a selectable 
mode of the Digitec console. 
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FIGURE 2 



Load Cells and Connections 



FIGURE 3 

Balance Control Box 
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III. SYSTEM KINEMATICS 



A. BASIC BALANCE STRUCTURE 

The basic balance structure consists of a frame that en- 
closes the balance mechanism and the primary support struc- 
ture upon which the model is mounted. 

The balance frame is a triangular base and two walled 
sides to which the internal balance mechanism is mounted. 

The base has castering wheels at each corner in order to 
provide a method of levelling the base of the balance. The 
top of the base serves as the floor of the balance. Two 
levels are attached to the floor and oriented at right angles 
to each other; the levels are used in conjunction with the 
levelling lugs in order to provide a horizontal- reference 
platform. The walls of the balance frame are inelastic 
under the maximum loads specified for the balance. The 
outside of the back wall serves as a mounting point for an 
electrical multiple outlet and a vibrator motor. The inside 
of the back wall supports various lever axes and flexure 
rod anchor points. 

The primary support member consists of a main column, 
a cross bar and model strut supports. The main column is a 
vertical hollow tube within which is mounted a vertical flex- 
ure for the transmission of lift forces. The main column 
serves as a method of transmitting model forces to the bal- 
ance. It also serves as a foundation for a convenient mod- 
el support system in the form of a cross bar to which the 
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span support struts are attached. A tail support arm is 
fastened to the upper area of the main column. This arm 
holds the tail support strut of the model, thus transmitting 
moment forces to the main column. An angle of attack motor 
can vary the position of the support arm so as to vary the 
angle of attack of the model. All support struts can be 
fastened to the cross bar and tail support arm so as to 
accommodate models of various proportions. 




1 1 
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PRIMARY SUPPORT COLUMN 



B. INTERNAL BALANCE MECHANISM 



The general principle of the MK I balance is measurement 
of the longitudinal forces required to restrain the movement 
of the primary support member when it is subjected to a 
loading condition through its attachment to a test model. 

The internal balance mechanism serves to restrain the motion 
of the primary support member in the lateral and longitudinal 
planes. These restraints are provided by flexure rods 
attached to the main column is such a way as to fix its 
position relative to the frame. The longitudinal and ver- 
tical flexure rods transmit their restraining forces through 
linkages to force flexures that create deflections in the 
respective load cells. The lateral flexure rods have no 
transmission capability. The linkage and flexure arrange- 
ment serves to reduce the loads applied to the load cells 
and to confine the internal mechanism to reasonable dimensions. 

The internal mechanism possesses a capability for some 
degree of adjustment in order to align the primary support 
member vertically. This adjustment is made possible by the 
use of threaded ends on both the flexure rods and force 
flexures. These adjustable ends provide variations in 
flexure length. 

All flexures are constructed so that their ends consist 
of a web section that assures the forces are transmitted 
axially along the flexure. The webs are susceptible to 
damage from loads that have shear or torsional components. 

It is therefore very important that longitudinal flexures 

i \ 
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INTERNAL MECHANISM SCHEMATIC 
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are oriented in a plane parallel to the back wall of the 
balance frame, and that the lateral flexures lie in a plane 
parallel to the front wall of the frame. All horizontal 
flexures must be parallel to the frame base. The anchor 
blocks to which one end of each force flexure is attached 
allow adjustment of the end in a plane perpendicular to the 
flexure axis such that the force flexures can be placed in 
the proper orientation. The remaining end of each force 
flexure and both ends of the restraint flexures are held in 
the proper position by their attachment to invariant fix- 
tures. These degrees of freedom are necessary in order to 
provide the capability of disassembly and assembly of the bal- 
ance as desired. When the above orientations are satisfied, 
all flexures are assured to carry their respective loads 
with no danger of damage to their fragile web ends. 

A complete alignment of the balance consists of a hori- 
zontal base, a vertical primary support column and the 
flexure orientation discussed in the previous paragraph. 

C. DERIVATION OF THE OUPUT MATRIX 

By evaluating the condition of the static equilibrium 
as applied to the primary support member, the forces trans- 
mitted to the three load cells can be derived in a simple 
matrix format that yields the loading condition imposed on 
the balance. The addition of an independent equation with 
angle of attack as a variable produces a four element vector 
that describes the longitudinal state of the model. 
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Equating static forces: 



L 



F 1 + F 2 + D " 0 
F 3 . L « 0 

Moments about strut pinnacle: 
M - F^a - 2 F 2 a = 0 
(a = 30.75 inches) 

Solving for restraint forces: 
F 1 = -(M/a + 2D) 

F 2 = (M/a + D) 

F 3 = - (L) 



Applying leverage relations 
to restraint forces produces 
the forces measured by the 



load cells: 

b/c = 2/3 d/d =1 f/e = 1/10 

P 1 = -2/3F 1 = 2/3 (M/a + 2D) 

P 2 = ' F 2 = '( M/a + D ) 



P 3 = -1/10F 3 = -1/10 (L) 













Collecting the above expressions 
into a convenient matrix format: 



1 


2/3 


4/3 


0 




N 

M/a 




M 




-2a/3 


4a/3 


0 




,p i 


= 


-1 


-1 


0 


< 


D 


> < 


D 


> = 


3/2 


1 


0 


< 


P 2\ 




•0 


0 


- 1/10 




i L . 




L 

J 




0 


0 


-10 




P 3 

V i 7 



Where fp} = [c]fF] [fJ = [t] fp} [t] = [cf 1 
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Since the state vector is in terms of the voltage outputs 

of the individual data channels, a conversion matrix must be 

formulated in order to produce an output vector in a more 

convenient format: 

/ \ 




M - (in-lbs) 
< D - (lbs) 

L - (lbs) 



v / 

by considering the slopes of the load cell outputs obtained 
from the calibration runs on the balance, a matrix was 
constructed that transformed load cell electrical signals 
into moment, drag and lift values in conventional units of 
measure. The method of obtaining the slopes of the cali- 
bration curves is explained in detail in Appendix B. 

For the horizontal force flexures: 



Pi = (A n ) M + (A 12 ) D 
P 2 - (A 21 ) m + (A 22 ) d 
F or the vertical force flexure: 

P 3 - (Ajj) L 

From the drag calibration curves: 

P x (mv) = +29.861 (D lbs) 
( mv ) = -22.656 (D lbs) 
From the moment calibration curves: 
P^ (mv) = .491 (M in-lbs) 
P 2 (mv) = -.730 (M in-lbs) 
From the lift calibration curves: 

P 3 = 4.25 (L lbs) 



23 



V 



Using this information, a coefficient matrix can be 



constructed : 



/ . A 

Ov) 




.491 29.861 0 




' M (in-lbs) ' 


/ P 2 (mv) 


>- 


-.730 -22.656 0 
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D (lbs) 


f — N 
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0 0 4.25 




L (lbs) 

/ 



M ■ DO- 1 M • H is the coefficient 

matrix 

(The transformation or output matrix.) 



[P} = [c]fF] 
[C]' 1 - [T] 



In conclusion, 
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IV. TESTING AND CALIBRATION PROCEDURE 



A. ALIGNMENT OF INTERNAL MECHANISM 

The alignment of all flexures is preferably carried out 
with the aid of a surveyor's transit and machinist's steel 
rule. The internal alignment should begin by blocking the 
wheels so that the balance will not accidentally roll, and 
then carefully levelling the base. The transit should first 
be placed so that the internal mechanism is viewed at ap- 
proximately a right angle to the back wall of the frame. 
After levelling the transit, the lens can be swung in azi- 
muth and elevation such that the horizontal grid of the 
transit passes through both webs of each horizontal force 
flexure. If this is not the case, then the flexure in 
question is not horizontal and the proper adjustment to the 
moveable end must be made. The vertical grid of the transit 
is used for alignment of the lift force flexure as explained 
in IV B. 

The steel rule should be used to measure the distance 
of the web ends of the two drag force flexures from the 
back wall of the balance. If a force flexure is not paral- 
lel to the back wall, its orientation must be corrected by 
adjusting its moveable end. This method assures that the 
drag force flexures, horizontal restraint flexures, lift 
force flexure and lift restraint flexure all lie in the 
longitudinal plane of the balance. 

With the transit in the same position relative to the 
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FIGURE 8 




Lever arms and anchor blocks 



FIGURE 9 
Load cells and 
Horizontal force 
flexures 




the main column. Adjustments to the primary support align- 
ment are made by varying the length of the two pairs of hori- 
zontal restraint flexures attached to the main column. It 
is imperative that the primary support be examined from 
planes that contain the restraint flexures, so that any 
adjustment to the column in one plane does not move the 
column in the other. 

When the above described alignment procedure has been 
completed, no further adjustment is necessary. As long as 
the balance base is levelled, the primary support member is 
vertical to the earth and therefore outputs will be accurate 
readings of the loads the balance is subjected to. 

C. TESTING OF ELECTRICAL CIRCUITRY 

All electrical apparatus for operation of the balance 
was connected to a common ground. The most readily available 
junction was the receptacle mounted to the balance back wall. 
The control box, vibrator, data acquisition system and angle 
of attack motor were joined to this receptacle. The recep- 
tacle in turn was connected to any available grounded elec- 
trical outlet. This method provided the desired common 
ground condition. 

The power/signal cables from the control box were attached 
to the appropriate transducers by means of quick disconnect 
fittings. The multiplex cable was utilized between the con- 
trol box and the data acquisition unit. 

Proper voltage sensing was checked by means of hand 
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applied forces at the center of the crossbar. By pushing 
down on the crossbar, the lift load cell was to produce a 
negative voltage. When a basically horizontal force was 
applied toward the rear of the balance, approximately a 
positive drag force, the top drag load cell was to produce 
a positive voltage while the bottom cell produced a nega- 
tive voltage. The angle of attack transducer was to pro- 
duce an increasing positive voltage for an increasing model 
pitch angle. When any applied force was removed from the 
balance, the load cell readings returned to their original 
values. This insured that there was no electrical drift 
from the bridges' zero readings, which could be set by use 
of the potentiometers mounted on the face of the control 
box . 



D. CALIBRATION OF THE BALANCE 
1 . General 

The balance had to be calibrated for loading con- 
ditions of lift, drag, pitching moment and combinations of 
these three independent conditions. Calibration for angle 
of attack was also required. The various loadings required 
were achieved through the use of pan weights or hydrau- 
lically applied forces. 

Loading the balance was accomplished through the 
use of a calibration bar designed for the purpose. The bar 
attached to the span support and tail struts much as a full 
model. The calibration bar had holes drilled into it such 
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that swage fittings could easily be attached to it. The swage 
fittings were mounted on steel cables that supported pan 
weights or hydraulic loads as appropriate. The cables were 
used because they were capable of being run over a pulley 
system in order to orient the applied loads horizontally 
or vertically as required. 

For each type of applied load, the following initial 
procedures were completed prior to the calibration. First, 
the balance was secured such that it would not move as a 
load was applied to it, then the balance was carefully lev- 
elled. Once this was done, the installed calibration bar 
was levelled using the angle of attack motor. This was most 
easily accomplished with a precision level that rested on 
the calibration bar. The appropriate loading devices were 
next attached to the calibration bar, the vibrator motor 
turned on, and all four output channels were brought to zero 
voltage readings with their respective potentiometers. 

The data for each calibration run were acquired in 
the following manner. First, a suitable incremental load 
was applied to the calibration bar. Ten-pound increments 
were used with pan weights and approximately fifty-pound 
increments when using the hydraulic method. If the load 
was by use of pan weights, the pan had to be steadied from 
its pendulum- like oscillations in order to obtain steady 
output readings. The Digitec multiplexing unit was then 
stepped through the three channels in order to read the 
digital voltmeter outputs and print out the values on the 
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FIGURE 10 



Balance Anchoring Procedure 



FIGURE 11 
Calibration bar 
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automatic recording device. Once the maximum desired load 
was reached, readings were also taken while reducing the 
load in increments in order to examine the data for evidence 
of hysteresis. The swage fittings were frequently checked 
for secure attachment to the calibration bar in order to 
prevent a serious accident due to the weights falling from 
their suspended position. 

2 . Lift Calibration 

Negative lift was applied through pan weights; the 
swage fitting of the pan was suspended from the lower center 
hole of the calibration bar so that the center of gravity 
of the load was applied along the lift restraint flexure. 

One hundred fifty pounds was the maximum practical load that 
could be applied in negative lift due to balance size con- 
siderations . 

Positive lift was applied hydraulically. The top 
center hole of the calibration bar was centered under the 
suspension point of the hydraulic piston used to create the 
loads. A plumb bob was required in order to assure exact 
alignment. A five’ hundred pound load cell was placed in the 
suspension system between the hydraulic piston and the cali- 
bration bar. The five hundred pound load cell was used to 
measure the load applied to the balance. Hydraulic pressure 
was supplied by a portable hand pump attached to the hydraulic 
piston. A small preload was first applied in order to bring 
the suspension system taut, then all balance output channels 
and the five hundred pound load cell bridge with its separate 
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DVM were zeroed before commencing the calibration. The 
maximum positive lift load applied was five hundred pounds , 
which transmitted the design maximum load of fifty pounds 
to the lift load cell. 

3 . Drag Calibration 

Drag loading was accomplished using pan weights up 
to a total of seventy-five pounds. A frame structure and 
pulley assembly was used to orient a cable such that the 
pan weights would apply a horizontal load to the balance. 

The swage fitting on the cable was pinned to the endmost 
hole on the calibration bar. The balance was positioned 
such that when the cable was suspended over the pulley, it 
was colinear with the longitudinal axis of the calibration 
bar to which it was attached, and therefore the drag force 
was being applied along the longitudinal axis of the balance. 
The pulley was then adjusted in elevation such that the 
cable was horizontal. The surveyor's transit was used to 
make this adjustment. 

4 . Moment Calibration 

The moment calibration procedure was carried out 
with the use of pan weights and hydraulic force. The bal- 
ance was positioned such that either of the holes that de- 
fined the upper four inch moment arms of the calibration 
bar were directly beneath the suspension point of the hy- 
draulic cylinder. The general procedure for creating and 
measuring the hydraulic force was the same as that used for 
lift calibration. 
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All moments up to four hundred inch pounds were 



applied by suspending pan weights from the appropriate 
moment arms of the calibration bar. Larger moments, up 
to the design maximum load of seventy-five foot pounds, had 
to be created hydraulically because suspending any more pan 
weights became physically impractical. Both the four inch 
and two inch moment arms were used to produce moment loadings 
when using suspended weights. This provided a cross check 
of balance alignment as discussed in Section VI. 

5 . Angle of Attack Calibration 

Calibration of the angle of attack transducer did 
not require that the balance be positioned or levelled in 
any particular manner. The pitch drive system was posi- 
tioned by means of its associated electrical controls at- 
tached to the forward wall of the balance. The controls 
encompassed start, stop, forward, reverse and speed func- 
tions. The geometric angle of the tail strut support arm 
was measured by a permanently installed rotary gauge on the 
back wall of the balance. Adjustable limit switches on 
this gauge serve to prevent damage to any particular model 
of exaggerated proportions due to contact with the balance 
as the model is rotated. 

The pitch drive was first run to its full nose 
down attitude, of -36 degrees, and the transducer output 
voltage was set to a minimum output reference. The pitch 
mechanism was then rotated in successive nose up increments 
of five degrees and the output voltage recorded; this was 
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continued until the maximum positive pitch angle of +36 






degrees was achieved. 
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FIGURE 14 







FIGURE 16 

Drag Calibration Procedure 
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FIGURE 17 

Drag Calibration Frame 
Closeup 
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FIGURE 18 




Negative moment 
4 inch arm 

Pan weight loading 



FIGURE 19 
Positive moment 
4 inch arm 
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FIGURE 20 






Combined loading 
procedure 




FIGURE 21 
Digitec data 
Acquisition console 
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Y, PRESENTATION OF DATA 



All calibration data are presented in tabular form. 
Each table of data is followed by a graphical plot except 
in the case of combined loadings. 



1 



\ 
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TABLE I 



LIFT CALIBRATION DATA 



LIFT LOAD 
pounds 




LOAD CELL OUTPUT 
millivolts 






P 1 


P 2 


P 3 


0 


0 


0 


0 


50 


0 


0 


212 


100 


0 


0 


423 


150 


0 


0 


635 


200 


0 


0 


855 


250 


1 


0 


1067 


300 


1 


0 


1276 


350 


2 


-1 


1489 


400 


3 


-1 


1702 


450 


3 


-2 


1915 


500 


3 


-2 


2125 


0 


0 


0 


0 


-12.5 


0 


0 


-52 


-25 


0 


0 


-105 


-37.5 


0 


0 


-158 


-50 


0 


0 


-212 


-60 


0 


0 


-254 


-70 


0 


0 


-297 


-80 


0 


0 


-339 


-90 


-1 


0 


-381 


-100 


-1 


0 


-424 


-125 


-1 


0 


-529 


-150 


-1 


0 


-636 
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LOAD CELL (P3) OUTPUT- ( millivolt a ) 




\ 
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LOAD CELL (P3) OUTPUT- (millivolts) 



TABLE II 



DRAG CALIBRATION DATA 



DRAG LOAD 
Pounds 




LOAD CELL OUTPUT 
Millivolts 






P, 


P, 


P 




1 


2 


3 


0 


0 


0 


0 


12.5 


375 


-282 


0 


25 


751 


-565 


0 


37.5 


1125 


-849 


0 


50 


1497 


-1132 


0 


60 


1792 


-1357 


0 


70 


' 2087 


-1582 


0 


75 


2235 


-1696 


-1 


70 


2098 


-1591 


0 


60 


1802 


-1365 


0 


50 


1508 


-1140 


0 


37.5 


1132 


-854 


0 


25 


756 


-571 


0 


12.5 


379 


-285 


0 


0 


0 


0 


0 
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POSITIVE DRAG CALIBRATION CURVE 



\ 
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LOAD CELLS (PI) a (P2) OUTPUT- (millivolts ) 



TABLE III 



MOMENT CALIBRATION DATA 





MOMENT 

Inch-Pounds 




LOAD CELL OUTPl 
Millivolts 


2 inch arm 


4 inch arm 


Pi 


P 2 


0 




-L 

0 


0 




0 


0 


0 


25 




12 


-18 


50 




24 


-36 




50 


24 


-36 


75 




36 


-55 


100 




48 


-73 




100 


49 


-73 


120 




58 


-88 




150 


74 


-110 




200 


98 


-147 




240 


118 


-177 




400 


196 


-292 




800 


391 


-586 




900 


445 


-657 




240 


118 


-176 




200 


98 


-147 




150 


74 


-110 


100 




48 


-73 




100 


49 


-73 


57 




36 


-55 


50 




24 


-36 




50 


25 


-36 


25 




12 


-18 


0 


- 


0 


0 




0 


0 


0 


0 




0 


0 




0 


0 


0 
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TABLE III (continued) 


2 inch arm 




4 inch arm 


P 1 


P 2 


-25 






-11 


18 


-50 






-24 


36 






-50 


-24 


36 


-75 






-36 


54 


-100 






-48 


73 






-100 


-49 


72 


-120 






-58 


87 






-150 


-74 


109 






-200 


-98 


146 






-240 


-119 


175 






-400 


-196 


291 






-800 


-391 


585 






-900 


-444 


656 






-240 


-119 


175 






-200 


-98 


146 






-150 


-74 


109 


-100 






-48 


73 






-100 


-49 


72 


-75 






-36 


54 


-50 






-24 


36 






-50 


-25 


36 


-25 






-11 


18 


0 






0 


0 






0 


0 


0 
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design limit- 75ft-lbs 
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LOAD CELLS (PI) S (P2) OUTPUT" ( millivolts ) 



0001 



O 

O 

00 




LOAD CELLS (PI)a(P2) OUTPUT- (millivolt*) 



TABLE IV 



ANGLE OF ATTACK CALIBRATION DATA 



ANGLE OF ATTACK 
Degrees 


TRANSDUCER OUTPUT 
Millivolts 


36° * 


P 4 

74 


35° 


126 


30° 


371 


25° 


590 


20° 


786 


15° 


984 


10° 


1160 


5° 


1335 


0° 


1493 


-5° 


1643 


-10° 


1784 


-15° 


1920 


-20° 


2050 


-25° 


2175 


1 

o 

o 


2290 


-35° 


2398 


o * 
-36° 


2420 



*Note: +36° is the physical limit of the AOA rotary rheostat; 

larger angles will ruin the transducer. 



< 
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ANGLE OF ATTACK - (degreet) 



FIGURE 27 

ANGLE OF ATTACK CALIBRATION CURVE 
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TABLE V 



COMBINED LOADING DATA 



LIFT LOAD 
Pounds 


DRAG LOAD 
Pounds 


MOMENT 

In-lbs 


LOAD 


CELL 

Millivolts 


OUTPUT 








P, 


Po 


P, 








1 


2 


3 


0 


0 


0 


0 


0 


0 


-12.5 


10 


-50 


289 


-207 


-52 


-25 


20 


-100 


578 


-417 


-106 


-37.5 


30 


-150 


865 


-625 


-159 


-50 


40 


-200 


1152 


-833 


-213 


-37.5 


30 


-150 


868 


-627 


-159 


-25 


20 


-100 


582 


-419 


-106 


-12.5 


10 


-50 


291 


-209 


-52 


0 


0 


0 


0 


0 


0 



VI. CONCLUSIONS 



A. ALIGNMENT ACCURACY 

The balance alignment described in Section III provided 
an accurate orientation for the flexures of the internal 
mechanism. The proper alignment of the flexures was required 
primarily to transmit restraint forces along the flexure 
axes, thus preventing damage to the flexure webs. The verti- 
cal alignment of the primary support was necessary in order 
to accurately resolve an applied load into its independent 
components of lift, drag and pitching moment. This alignment 
was in fact very accurate due to the careful levelling of the 
balance and transit as well as the use of significant colinear 
points to ascertain the necessary orientations. Once the 
balance was aligned by use of the transit, the most readily 
available crosscheck as to the resolution capability of the 
balance was to apply independent loadings of moment, lift 
and drag and observe data channels for indications of coupled 
output . 

If the primary support was vertical, then a pure lift 
load applied to the balance would create an output only on 
the lift data channel. Any output on the drag channels 
would be evidence that the column was not vertical and conse- 
quently the vertically applied load was being resolved such 
that it produced a drag component directly proportional to 
the angular misalignment of the primary support, as well as 
creating a moment about the bottom anchor point of the lift 



54 



restraint flexure. By the same reasoning, an isolated 
drag load would produce some output on the lift channel 
that was proportional to that angular misalignment. Posi- 
tive and negative moments of the same absolute value should 
produce output voltages of the same absolute magnitude. An 
additional method of cross checking the primary support 
alignment by means of applied pitching moments consisted of 
examining equal moments created by suspension from different 
moment arms. If the center of the calibration bar did not 
lie directly above the axis of the lift restraint flexure, 
the moment output channels would show different values for 
a given moment depending on whether the weights were sus- 
pended from a two or four inch moment arm. Because of the 
above considerations, an initial moment calibration should 
be made after the balance is aligned. This is a relatively 
simple process since it requires no auxiliary framework to 
apply the loads. Any minor adjustment that is required in 
order to produce consistent moment output signals must be 
made by changing the lengths of the two horizontal restraint 
flexures utilizing their adjustable ends. The top restraint 
flexure will change the angle that the lift restraint flexure 
makes with the vertical, since the bottom of the vertical 
flexure is anchored outside the main column. The bottom 
horizontal restraint flexure will change the horizontal 
location of the mounting eyes of the span supports. Once 
the initial moment calibration is satisfactory, the complete 
procedure can be carried out with a high degree of confidence 
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that the balance is correctly aligned. 

B. DATA ANALYSIS 

The graphical presentation of calibration data contained 
in Section V illustrates the fact that the load cell output 
is linear within the design load ranges of the balance. A 
least squares fit for all applicable data was conducted by 
use of a computer program which is described in detail in 
Appendix B. The results of this program confirm the graphi- 
cal evidence of balance linearity. The only exception to 
this is the angle of attack transducer; however, since this 
variable is an independent factor, it has no effect on force 
resuolution by the balance. In addition, a least squares 
second order polynomial approximation of the transducer 
output produced a very accurate representation of the angle 
of attack calibration curve. The maximum error between the 
two curves was less than one percent: 

a a a 

P 4 CAL DATA LEAST SQUARES ERROR 

(1160 mv) (10.0 deg.) (10.09 deg.) (.0925 deg.) 

The vibrator motor made very little difference in output 
signals. This indicated the virtual absence of mechanical 
hysteresis in the balance. The calibration of the balance 
does utilize the data acquired with the vibrator on, because 
in each case it was slightly more exact as applied to linear 
output . 

The calibration data for lift, drag and moment were also 



\ 
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analyzed for deviation from computed mean values. The 
results appear in the output of the data analysis computer 
program as standard deviations for each load cell under the 
three conditions of independent loading. A correlation of 
these deviations produces the following results: 

DRAG 









7.1 



3.55 



Utilizing the balance output matrix for determination of 
drag : 

^DRAG = (.0684)0?! + (-046)0^ = .649 lbs 
MOMENT 



R 



OMENT 



Opi 


f.865' 


fe P2 

V J 


r 621 , 



(^MOMENT (2.122 30^! + ( 2.79 7)0^2 -3.572 in lbs 



LIFT 



fell fep2 1 I 1 * 69 } 

l Jlift j 



0LIFT = (.235)0^3 = .397 lbs 
For typical test conditions, S = 3 ft' 



q = 20 lb/ft ; therefore, 



CTc d - 



tfc M ■ 



(.649 lbs) 

(20 lb/ft 2 ) (3 ft 2 ) 



= .0108 



(-3.572 in-lb) 

T 



T 



c = .75 f t , 



.0066 



cr. 



C L * 



(20 lb/ft )(. 75 ft) (3 ft ) (12 in/ft) 
(.397 lbs) 



(20 lb/ft 2 ) (3 ft 2 ) 



= .0066 
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The output matrix calculated a loading vector that had 
the values of the independent loads computed accurately. 
However, unacceptable values for pitching moment were pro- 
duced from an independent drag load as well as slight amounts 
of drag from a pure moment application. These aberrations 
are displayed in the output section of the data analysis 
computer program. They were a direct result of the trans- 
formation matrix coupling all deviations in horizontal load 
cell signals and apportioning the deviations erroneously. 

The primary cause of such distribution is the fact that the 
voltage full scale range for drag was an order of magnitude 
greater than that for pitching moment. Since the balance- 
computed loading vector is the product of a linear system, 
the computational errors committed through the use of the 
output matrix can be analyzed by the principle of superposition 



(f + AF] = [t] fp + AP} (f) = [t] [ap} 



'am' 




-2.122 


-2.797 


0 




A p i' 


Jad 


> = 


.0684 


.046 


0 


< 


AP 2 l 


AL 

V J 




0 


0 


. 235 




AP 3 



AM (in- lbs) = -2.122AP 1 - 2.797 A P £ 

AD (lbs) = .0684AP 1 + .046AP 2 

It is clear that small errors in signals P^ and P 2 will 
cause unacceptable error probabilities. A standard deviation 
in drag coefficient of .0108 is of the same order as the 
profile drag of a test model. Obviously, a standard deviation 
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in pitching moment of over three inch-pounds is excessive 
for the purposes of model testing. In its present con- 
figuration, the balance could not resolve combined loading 
conditions of drag and moment; the drag resolution was 
reasonable but the moment measure was not: 

COMBINED LOAD OUTPUT BALANCE RESOLUTION 

M = -200 in-lbs P x = 1152 mv M = -114.7 in-lbs 

D = +40 lbs P 2 = -833 mv D = +40.46 lbs 

C. RECOMMENDATIONS 

As stated previously, the primary reason for the failure 
of the balance to resolve combined loadings of moment and 
drag is the scale range differences in load cell output be- 
tween drag and moment loadings. The cause of the scale 
range incompatibility is the length of the span support struts 
? 1 + F 2 + D = 0 
-F^a - F 2 (a + b) + M = 0 

solving for restraint forces: 

F 1 = [£» + b ) - m] 

F 7 = 1 (Da + M) 

1 b 

If length a is decreased the re- 
lative magnitude of the moment 
and drag contributions to the 
restraint forces become more 
compatible. However, the prac- 
ticality of such a modification is limited by the necessity 
of having struts that are long enough to suspend a model in a 
tunnel test section. 
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APPENDIX A 



LOAD CELL DESIGN AND SPECIFICATION 



GENERAL 



After considering the geometric restrictions of the 
linkages comprising the internal balance mechanism as well 
as the desired load measurement accuracy, the primary 
specifications for the load cells were established. These 
parameters were a maximum dimension of three inches along 
the major axis of each cell and a maximum allowable deflection 
range of plus or minus eight thousandths of one inch. All 
remaining dimensions except the web thickness were derived 
from similar practical considerations. 'Aluminum was chosen 
for compatibility with existing balance construction. 

STATIC ANALYSIS 

Consider a load cell as an indeterminant structure com- 
posed of two sets of identical components as shown in figure 



M(x) = (M - Px); 

By the theorem of Castigliano, the elastic strain energy 



Since 1^, the total elastic energy stored can be 

approximated by 



( 29 ): 



The moment function in Piece 1 can be expressed as 
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TABLE VI 



LOAD CELL DIMENSIONS 



DIMENSION 


50 LB LOAD CELL 


140 LB LOAD CELLS 


a 


• 

o 

o 


• 

O 

o 


b 


1.882" 


1.98" 


c 


2.280" 


2.280" 


d 


1.140" 


1.140" 


e 


1.640" 


1.640" 


f 


0.121" 


0.170" 


g 


2.522" 


2.620" 


h 


0.50" 


0.50" 


j 


0.495" 


0.495" 


k 


10-32 clearance 


10-32 clearance 
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LOAD CELL STATICS 



FIGURE 29 



(1) cr-M 

J o • 

jv =2s, r*.°; - a to- [/>-■£ (s-;)J 

The moment function for component (2) is a constant value, 
M; 

^irfa ^ ft-) ■-$ -jfa -Al 

^ © 

for compatibility, Q, ~ ~ Q 2 » therefore 

( 2 ) 

#*£/ 0-4/4) =r 

3 ^ h/ ^ufaz) sU+^0-^Uj\ 



The stress on the surfaces of component (2) can be 
calculated : 



■ si vm)£_. /St ) 

I, ’ P '■r/ax, 



For aluminum, 



( 3 ) 




£-6 =■ /£?/£ s sLx 
e 7 £ 







APPLICATION 

The difference between the fifty pound load cell and 
the one hundred pound load cells is reduced to a required 
variation in web thickness in order to limit the 
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resultant deflection due to maximum applied load. 

For one hundred forty lb cell: 

Deflection of component (T) at point of load application 



is 



0 s AU-J) 

M> =^J J 



dP 



- 

2EL 






The total deflection of the cell along the load application 
axis is twice this value 




To calculate the minimum web thickness t: 
T- UJl! 7 /AT 



/a 






PJ t[/ rZO-M)] 



using the approximated dimensions 

/ = 2.4 inches h = 1.7 inches 

d = . 50 inches ^4^,= 16 x 10 ^ inches at 140 lbs load 

w = .495 inches 

&UL = . 5 = .20833 (I -dll) = .791666 /(| -d //) = 1.9 

2.4 



65 



i 



"d |f^ 



P (1 - d//) 3 = 6.859 h + /(I - d//) = 10.459 

t 3 = 140 6.859 (12.5) = 

(4) (.495) CIO 7 ) (16 x 10' 3 ) 10.459 

(4.4192 x 10" 4 ) (8.19748) = 3.62263 x 10‘ 3 ) 






.154 inches 



For an actual web thickness of .170 inches, the use of 
equation (3) gives the output voltage rates for the load 
cell connected to a balanced bridge circuit: 



£ = £ 5 t (1 - d//) 2 



I = wt 3 = (.495) (3.62263) x 10 -3 



1601 [h+/(l-d/,£)] 



12 



1.494334 x 10' 4 in 4 



e = 


(2 . 62) 3 


(.170) x 10 4 " 




.6547 


p 


(160) (1. 


4943) 




(1.98 + 2.119) 

v j 



12 



= 20 . 42 ^in/ in = 
lb 



2859 mv at 140 lbs 



The strain gauges were mounted on the side components 2 
of the load cells because the moment along these members does 
not vary with strain gauge positioning. 



For an excitation voltage of 5 volts DC and a gauge 
factor of 2.0, 1000y^ in/in. = 10 MV output from a balanced 
bridge circuit. 

x P =20.42 in/in x 140 lbs = 2859 MV maximum range; 
max 

lb 

this agrees closely with calibration values. 
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APPENDIX B 



DATA ANALYSIS COMPUTER PROGRAM 



GENERAL 

All data analysis was carried out through a FORTRAN IV 
program utilizing the Naval Postgraduate School IBM 360/67 
digital computer. The primary impetus for computer analysis 
resulted from the necessity to obtain accurate slope values 
from the various load cell calibration curves. Although 
reasonable results could be obtained by graphical methods, 
more precise information was available through the use of a 
least squares fit. 

An initial graphical plot of calibration data showed 
the balance output to be linear with respect to independent 
loading conditions; consequently, the least squares fit of 
data for each load cell was assumed to conform to a first 
order polynomial of the form y = ax + b where (a) represents 
the slope of the particular calibration curve and (b) the 
intercept of the curve along the ordinate. 

Calibration data for the angle of attack transducer 
necessitated the use of a higher order polynomial. The most 
functional ..solution to this problem was to use a packaged 
subroutine to calculate the coefficients of a fourth order 
polynomial, since a trial procedure indicated a good fit for 
such an order. It was estimated that an acceptable accurate 
expression for the model angle of attack could be produced 
by this method. 
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PROGRAM DETAILS 



All calibration data was arranged to be read in as a 
matrix INPUT (102,5) arranged in an easily accessible 
pattern as follows: 



l 1 2 3 4 5 

52 
67 
85 
102 

Two original subroutines C0MP1 and COMP2 were used to compute 
least squares first order polynomials for applicable load 
cells under independent conditions of moment, drag or lift 
loadings. C0MP2 is used only for analyzing lift data and is 
a necessity in order to avoid cumbersome linkage problems 
with the main program. The resultant slopes and intercepts 
are stored in matrices COEFF (3,3) and INCPT (3,3) for 
further processing. The format and information contained 
in these matrices is readily discernible from the form of 
the offline printer output. 

Once these steps have been completed, the main program 
produces the balance output matrix TRANS (3,3). The develop- 
ment and purpose of this matrix was discussed in detail in 
Section III. TRANS (3,3) is utilized by a third original 



MOMENT 

LOADS 


P 1 


P 2 


NULL 


NULL 


DRAG 

LOADS 


p i 


P 2 


NULL 


NULL 


LIFT 

LOADS 


NULL 


P 3 


NULL 


AOA 


NULL 


NULL 


P 4 
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subroutine, C0MP3. 

C0MP3 performs three major operations on the calibration 
data for moment, drag and lift. First, COEFF (3,3) is used 
to calculate the mean values of load cell output in milli- 
volts for each value of independent variable available. 

The second operation performed is the subtraction of these 
mean values from the actual calibration values in order to 
obtain a standard deviation for each load cell under each 
applicable loading condition. The standard deviations, 
expressed in units of millivolts, are stored in matrix 
SIGMA (3,3) for output purposes. The final operation of 
C0MP3 is to multiple the actual calibration values of load 
cell output by the output matrix in order to obtain the 
calculated loading conditions that would be produced by 
the balance. These computed loadings are stored as actual 
values in matrix ACTVL. 

A prepackaged subrouting LSQPL2 was used to evaluate 
the coefficients of the fourth order polynomial that expresses 
angle of attack as a function of transducer voltage. Some 
of the printed output from this subroutine was supressed 
for brevity and clarity. 

The printed output from the data analysis program 
arranges all calculated and indexing information in con- 
venient reference form. The program listing contains adequate 
comment to enable the interested reader to follow the sequence 
of events contained in the program. 
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MOMENT DATA ANALYSIS 
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COMPARATIVE OUTPUTS COMPUTED LOAD 
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ANGLE OF ATTACK DATA ANALYSIS 
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SLCFES INTERCEPTS STANDARD DEVIATIONS 
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on ooooooo Hoooonoooo ui ooooonooo oooonoooooon 



ccccccccccccccccccccccccccccccc 



c 

c 



c 

c 

c 

c 

c 



C DATA REDUCTION PROGRAM 



c 

c 



15 



ccccccccccccccccccccccccccccccc 



INTEGERS A t G 

REAL*4 MEAN , INPUT ,INCPT ,ACTVL 
REAL*8 X,F2,WI,Y,DELV,B, SB, TITLE 
DIMENSION X(17),F2117),*I(17),Y(17),CELY(17), 
XE (11 ) ,S8( 11 ), C0£FF(3, 2 ), TRANS (3,3) , SIG(3,3) , 

X I NC P T (3,3) , I N PUT ( 102 » 5 ) , POINTS (60,5) , 

XMEAN( 102,3) ,ACTVL( 102,2) 

COMMON INPUT 



CCCCCCCCCCCCCCCCCC 
C C 
C READ IN DATA C 
C C 
CCCCCCCCCCCCCCCCCC 



DC 5 1=1,102 

READ (5,500) ( I NPUT ( I , J ) , J = 1 , 5 ) 

CONTINUE 

CO 35 A= 1 , 3 

GO TO (10,15,20) » A 



CCCCCCCCCCCCCCCCCCCCCCCCC 
C C 
C MOMENT CALCULATIONS C 
C C 
CCCCCCCCCCCCCCCCCCCCCCCCC 



K = 52 
L = 2 
M = C 

GC TO 25 



CCCCCCCCCCCCCCCCCCCCCCC 
C C 
C DRAG CALCULATIONS C 
C c 
CCCCCCCCCCCCCCCCCCCCCCC 



K = 1 5 
L = 2 
M = 52 

GC TO 25 
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c 

c 

c 

c 
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20 
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c 
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c 
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c 
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25 
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C 
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c 

c 

c 

c 

c 

20 

35 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



ccccccccccccc ccccccccccc 

C LIFT CALCULATIONS C 
CCCCCCCCCCCCC CCCCCCCCCCC 



M = 6 7 

GC TC 30 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



c c 

C CCNP1 PERFORMS A LEAST SQUARES OPERATION ON C 
C TEE TEST DATA FOR MOMENT AND DRAG; IT THEN C 
C FRCVICES THE SLOPES ANC INTERCEPTS OF THE C 
C ASSUMED STRAIGHT LINE RELATION BETWEEN C 
C APPLIED LOAC AND LOAD CELL OUTPUT. C 
C C 



CCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCC CCCCCCCCCCCCCC 



CALL CCMPK K,L,M, A,COEFF, I NCP T , POINTS ) 
GC TC 35 



CCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCC CCCCCCCCCCC 



c c 

C CCNP2 PERFORMS THE IDENTICAL OPERATIONS C 

C CAFRIEC OUT BY COMPlJ HOWEVER, IT OPERATES C 

C ONLY ON LIFT DATA IN ORDER TO SIMPLIFY THE C 
C OVERALL OPERATIONS OF CCMP1. C 

c c 



ccccccccccc ccccccccccc-: ccccccccccccccccccccccccc 



CALL CCMP2( K, L, M, COEFF, INCPT, PO INTS ) 
CONTINUE 



CCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCC CCCCCCCCCCC 



c c 

C TRANS ( I » J) IS THE BALANCE OUTPUT MATRIX; IT C 
C CONVERTS LOAD CELL OUTPUT IN MILLIVOLTS TO C 
C MCMENT, DRAG AND LIFT. C 
C C 



CCCCCCCCCCCCC CCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCC 



CET=CCEFF(1,1 ) *CQ EF F ( 2,2)-C0EFF( 2, 1)*C0EFF( 1,2) 

TRANS (1,1) =CCEFF(2,2 )/DET 

TRANS ( 2,2 ) =CQEFF ( 1, 1J/DET 

TRANS (1,2) =-CCEFF( 1,2)/ GET 

TRANS(2,1) =-CCEFF(2»l)/0ET 

TRANS (3,1 )=C.O 

TRANS(3,2)=0.0 

TRANS(3»3) =CCEFF(3»3 ) 

DO 60 G = 1 » 3 

GC TC (40,45,50) , G 



77 



i 



ccccccccccccccccccccccccc 
c c 

C MOMENT CALCULATIONS C 
C C 
CCCCCCCCCCCCCCCCCCCCCCCCC 



L = 2 
M ” C 

WRITE (6,501 ) 
WRITE (6,502) 
WRITE (6,503) 
WRITE (6,504) 
60 TO 55 



CCCCCCCCCCCCCCCCCCCCCCC 
C C 

C DRAG CALCULATIONS C 

c c 

CCCCCCCCCCCCCCCCCCCCCCC 



K = 15 
L = 2 

|Yj — G 2 

WRITE (6,505) 
WRITE (6,502) 
WRITE (6,506) 
WRITE (6,507) 
GO TO 55 



CCCCCCCCCCCCCCCCCCCCCCCC 
C C 
C LIFT CALCULATIONS C 
C C 
CCCCCCCCCCCCCCCCCCCCCCCC 



K = 18 
L = 2 
M =67 

WRITE (6,508) 
WRITE (6,505) 
WRITE (6,510) 
WRITE (6,511) 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

C CCMP3 CALCULATES MEAN LCAC CELL OUTPUT 
C VALUES AND BALANCE LOADING VALUES BASED ON 
C THE SLOPES OBTAINED FROM COMP 1 AND CCMP2. 

C TFE SLB PROGRAM ALSO CCMPARtS MEAN OUTPUTS 
C WITH ACTUAL OUTPUTS IN ORDER TO CALCULATE A 
C STANDARD DEVIATION FOR EACH CELL UNDER 
C INDEPENDENT LOADINGS. 

C 

CCCCCCCCCCCCCCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCC 



OOOOoOOOO 



55 

60 



65 



C 

C 

c 

c 

c 

c 

c 

c 

c 

500 

501 

502 

502 

504 



505 

506 

507 



508 

50S 

510 

511 

512 
512 

514 

515 



516 

517 

518 



CALL CCMP3 (K ,L,M, COEFF, INCPT,S I G , ACT VL , T RAN S , 
XG, POINTS, DELTA, MEAN) 

CONTINUE 
WRITE (6,512) 

DO 65 1=1,17 
W I ( I ) =1 . 

F 2 ( I ) =INPUT (1 + 85,1) 

X ( I ) = INPUT ( 1 + 85,5) 

CONTINUE 
M = 17 
KM=-4 

CALL LSCPL2 ( M, KM,X, F2 , WI , Y, CELY , B, SB , T ITLE ) 
WRITE (6,512) 

WRITE (6,514) 

WRITE (6,515) 

WRITE (6,516) (COEFF (1,J) , J = 1 , 3 ) »(INCFT(1,J) , 
X J = 1 , 3 ) ,(SIG(1,J),J=1,3) 

WRITE (6,517) (COEFF (2 , J) , J=1 ,3 ) , ( INCPT (2, J ) , 
XJ-1,3) , (^> I G ( 2 , J ) , J = 1 • 3 ) 

WRITE (6,518) (CCEFF(3,J),J=1,3),(INCPT(3,J), 

X J = 1 , 3 ) , ( S I G ( 3 , J ) , J = 1 , 3 ) 



CCCCCCCCCCCCCC 
C C 

C FORMAT S C 

c c 

CCCCCCCCCCCCCC 



FGFMAT (5F10.2) 

FORMAT (• O' ,10 (/) ,T32 ,' MOMENT DATA ANALYSIS',///) 
FORMAT ( 'O' ,T32, 'COMPARATIVE OUTPUTS', 

XT74 , ' CCNPUT ED LOADING') 

FORMAT ( • O' ,///,T4, ' MOMENT • ,T24 ,'P1(MV)',T52, 
X'P2(MV)'»T72» 'MOMENT • ,T84, 'DRAG • ) 

FORMAT ('O' , / , T3 , ' (IN-LES)' ,T19, 'DATA ' , T 3 1 , 

X 'MEAN • ,T47, 'DATA' ,T59 , ' MEAN ' ,T7 3 , ' (IN-LBS) ' , 

XT85, ' (LBS ) ' ,// ) 

FORMAT ( ' O' , 10 ( / ) ,T33 , ' DRAG DATA ANALYSIS',///) 

FORMAT C0',///,T5,' DRAG' , T23 , ' PI ( M V ) • ,T52 , 
X'F2lMV)',T72,' MOMENT ' , T 84 , ' OR AG ' ) 

FORMAT ('O' ,/,T4,'(LSS)' , T1 9 , • DATA • , T3 1 , ' ME AN ' , 

XT 47, 'CATA' , T59 , ' MEAN ' , T7 3 , ' ( I N- LBS ) ' , 

XT 8 5 , ' (LBS) ' ,// ) 

FORMAT ( 'O' , 10 ( /) ,T28,* LIFT DATA ANALYSIS',///) 

FORMAT ('O' »T30, 'COMPARATIVE OUTPUTS', 

XT6C , ' COMPUTED LOAD* ) 

FORMAT ( ' O' ,///,T6, ' LIFT' ,T35»'P3(MV) ' »To6,' LIFT' ) 
FORMAT CO* ,/ ,T5, ' (LBS )• ,T29, 'DATA* ,T44, 'MEAN • , 

XT 6 5 , ' (LBS) ' ,//) 

FORMAT (• C' , 10 ( / ) ,130, 'ANGLE OF ATTACK DATA ANALYSIS' 
FORMAT ('O' , 5 ( / ) ,T36, • CATA FITTING VALUES',/) 

FORMAT (• O' ,T16, ' SLOPES' ,T45, ' INTERCEPTS' ,T70, 
X'STANlARD DEVIATIONS',//) 

FORMAT ( 'O' , 2 ( / ) , T 9 » ' M ' ,T 1 9 , ' D ' ,T29 , ' L ' , 
XT29,'M',T49,'D',T59,'L',T65,'M',T79,'C' ,T89,'L' ) 
FORMAT ('O' ,/ ,T2 , ' PI ' , 9F10 . 3 ) 

FORMAT (• O' ,/,T2,'P2' ,9F10.3) 

FORMAT ('O' ,T2,'P3',9F10.3) 

RETURN 

ENC 
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SLERGLTINE COMF1 ( K, L , M , A , CCEFF, INCPT, POINTS ) 
INTEGER** A 
REAL** INPUT, INCPT 

Cl PENSION P0INTS(K,L) , I NPUT (102 ,5 ) , CCEFF (3, 3) , 
X INCPT (3,3) 

CCMMCN INPUT 



cccccccccccccccc 
c c 

C INITIALIZE C 
C C 
CCCCCCCCCCCCCCCC 



SUMY = 0 . 
SLMP1 =0. 
SLMP2=0. 
SMSC1=0. 
SMSQ2=0. 
SM YP1-0. 
SMY P2 = 0 . 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

C FILL CATA SUBMATRIX PCINTS(I,J) C 
C C 

ccccccccccccccccccccccccccccccccccccc 



DC 110 1=1, K 
M = N + 1 

PC INTS (1,1 ) = I N PUT (M, 1 ) 
DC 105 J=2 » L 
PC INTS ( I , J) =INPUT (M, J ) 
05 CONTINUE 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

C CARRY OUT LEAST SQUARES OPERATIONS C 

c c 

cccccccc cccccccccccccccc cccccccccccccccc 



SUMY = SUMY+PCINTS( 1,1 ) 
SLMP1=SUMP1+P0INTS( I ,2 ) 
SUMP2=SUMP2+P0INTS( 1,3) 

SGPl = FCIi\TS (I ,2 ) **2 
SQP2= PO INTS ( I , 3 ) **2 
SMSG1=SMSC1+SCP1 
SM.SG2 = SMSC2 + SCP2 
YP1 = PCINTS( I,1)*P0INTS(I,2) 
YP2=PCINTS(I,1)*P0INTS(I,3) 
Si w YP1 = SMYP1 + YP1 
SMYP2=SMYP2+YP2 
1C CONTINUE 



CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C C 
C STORE SLOPES AND INTERCEPTS C 
C C 

ccccccccccccccccccccccccccccccccc 



IN CC EFF ( I , J ) AND INCPT(I,J) 

DNM1= (K*SMSCl-( SUMP 1**2) ) 
CNM2=(K*SMSC2-( SUMP 2**2) ) 

SLCPE1= (K*SMYP1-SUMP1*SUMY l/DNMl 
SLOPE 2= (K*SMYP2-SUMP2*SUMY) /DNM2 
CLT1= (SUMY*SM.SC1-SMYP1*SUMP1)/CNM1 
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CUT2= (SUMY*SMSC2-SMYP2*SUMP2»/DNM2 

CCEFF (1,A) =1/SLCPE1 

CGEFF(2,A)=1/SLGPE2 

CCEFF(3,A)=0.0 

I NCP T ( 1 ,A ) =CUT 1 

INCPT ( 2 , A ) = CUT 2 

INCPT(3,A) =0.0 

RETURN 

ENC 



SUBROUTINE CCMP2 ( K , L , M , C C E FF , I NC PT , POINTS ) 

R E AL*4 INPUT, INCPT 

CIMENSICN PCINTS ( K, L ) , INPUT (102, 5) ,CGEFF(3,3) , 
X I NCP T (3 ,3) 

COMMON INPUT 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C 

C SAFE OPERATIONAL STEPS AS SUBROUTINE CuMPl 
C 

CCCC COCO COCO CCCCCCCCCCCCCCCCCCCCCCCCC COCO CCCCCCC 



K 



SUMY=C. 

SUMP2=C. 

SMSG3=0. 
smyp3=o. 

CO 205 1=1 
M = M + 1 

PCINTS(I ,1) =INPUT(M,1) 

PO INTS ( 1 , 2) =INPUT (M, A ) 

SU MY = SUMY+ P C I NTS (1,1) 
SUMP2=SUMP3+P0INTS(I ,2) 

SUP3= PC INTS ( I , 2)**2 

SNSC3=SMSC3+SCP3 

YP2=PCINTS( I , 1)*PGINTS(I ,2) 

SMYP3=SMYP3+YP3 

CONTINUE 

CNM3= (K*SMSQ3-(SUMP3**2 ) ) 

SLCPE3= (K*SNYP3-SUMP3*SUMY J/DNM3 

C U73=( SUMY* SMS C3-SMYP3*SUMP3 ) /DNM3 

CCEFF( 1,3 )=0.0 

CCEFF(2,3)=0.0 

CDEFF(3,3) =SLCPE3 

INCPT ( 1 , 3 ) = 0 . 0 

INCPT(2,3) =0.0 

INCPT (3,3) =CUT 3 

RETURN 

END 



SUBROUTINE C0MP3 ( K , L , M , C CEFF , I NCPT , S I G , ACTV L , 
XT PANS, G,PC INTS , CELT A, MEAN ) 

REAL*4 INCPT, MEAN, INPUT 
INTEGERS G 

C IN ENSIGN ACTVL(K,3) , PC INTS (K , L ) , 0 EL T A ( K , 2 ) , 
XCOEFF (3,3) , INCPT (3,3) , SIG( 3,3) , MEAN (K, 3) , 
XTRANS (3,3) , INPoT< 102, 5) 

COMMON INPUT 
DELS M 1=0. 

DELSM 2 = 0. 



CCCC C CCCCCCCC CCCCCCC CCCC CCCC CCCCCCCCC 
C C 
C FILL CATA SUBMATRIX PCINTS(I,J) C 
C C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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205 CO 33 C 1 = 1, K 
N = N + 1 

PC IMS ( I ,1) =1 NPUT ( M» 1 ) 

I F ( G . EC .3 ) GO TO 310 
PCINTS(I,2)=INPUT(M,2 ) 

PCIN7S(I,3)=INPUT(M,3) 

GO TO 315 

21C PC1NTSU ,2) =INPUT(M,4) 

C 

C 

C CCCCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCC 

c c c 

C C CALCULATE MEAN LOAC CELL OUTPUTS C 

C C ANC BALANCE-COMPUTED LOADINGS C 

c c c 

c cccccccccccccccccccccccccccccccccccccc 

c 
c 

315 M E AN ( I,1) = PCINTS( I, 1) 

I F (G . EC .3 ) GC TC 320 

MEAN! I ,2)=PCINTS( I , 1 ) *CCEFF (1,G) 

MEAN< I ,3)=PCINTS( I, 1 )*CQEFF<2, G) 

ACTVL(I»1)=FCINTS(I»2)*TRANS(1,1)+P0INTS(I,3)*TFANS(1,2) 
ACTVL( I,2) = PGINTS(I , 2 IMRANS < 2, 1H-PCINTSU , 3 ) *T F AN S < 2 , 2 ) 
WRITE (6, 530) MEAN ( I , 1 ) , PCI NTS ( 1 , 2 ) , M£AN< 1 , 2 ) , PC INTS ( I , 3 ) , 
XMEAN (I ,3) , ( ACTVL( I , J) , J=l,2) 

530 FORMAT ( F8 . 1 , 3X , 2F 12 . 3 , 5X , 2F 12. 3, 5X , 2F 10. 3 ) 

GC TO 325 

320 M,EAN( I , Z) = PC I NTS ( I , 1 ) / COEFF < 3 , 3 ) 

ACTV L ( I , 1 ) = PC I NTS ( I, 2 ) *CC E FF ( 3 * 3) 

WRITE (6,531) MEAN( I , 1) , POINTS ( I ,2 ) ,MEAN( I, 2 ) , ACTVL (I , 1) 

531 FORMAT ( F10 .3 , 10X, 2F15.4 , 1CX, F10.3) 

32 5 CELT A ( 1,1 ) = POINTS (I, 2) -MEAN (I, 2) 

CELSG 1=DELT A ( I ,1)**2 
CELSM1=CELSM1+DELSQ1 
33C CONTINUE 

RG1=DELSM1/K 
SIG1=SCRT (RG1) 

GC TC ( 335,340,345) , G 
335 SIG(1»1)=SIG1 

GC TC 350 

34C S I G ( 1 ,2) =SI G1 

GO TO 350 

345 S I G ( 3 , 3 ) =S I G1 
C-0 TO 355 

350 DO 355 1 = 1 , K 

DELTA ( I ,2) =FCIN'TS(I ,3) -MEAN (I ,3 ) 

CELSQ2=DELTA< I, 2)**2 
DELSM2=CELSM2+DELSQ2 
255 CONTINUE 

P G2 = D ELSM2/ K 
S I G2 = S CRT ( R.G2 ) 

IF(G.EQ.l) S I G ( 2 , 1 ) = S I G 2 

IF(G.EC.2) SIG(2,2)=S IG2 

SIG(3,1)=0.0 

S I G ( 3 , 2 ) =0. C 

S I G ( 1,2) = 0.C 

SIG(2,3)=0.0 

RETURN 

ENC 
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